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ABSTRACT 

We use Sloan Digital Sky Survey (SDSS) data to investigate galaxy cluster prop- 
erties of systems first detected within DPOSS. With the high quality photometry of 
SDSS we derived new photometric redshifts and estimated richness and optical lumi- 
nosity. For a subset of low redshift {z ^ 0.1) clusters, we have used SDSS spectroscopic 
data to identify groups in redshift space in the region of each cluster, complemented 
with massive systems from the literature to assure the continuous mass sampling. A 
method to remove interlopers is applied, and a virial analysis is performed resulting 
in estimates of velocity dispersion, mass, and a physical radius for each low-z sys- 
tem. We discuss the choice of maximum radius and luminosity range in the dynamical 
analysis, showing that a spectroscopic survey must be complete to at least M* -f- 1 if 
one wishes to obtain accurate and unbiased estimates of velocity dispersion and mass. 
We have measured X-ray luminosity for all clusters using archival data from RASS. 
For a smaller subset (twenty-one clusters) we selected temperature measures from the 
literature and estimated mass from the M-Tj^ relation, finding that they show good 
agreement with the virial estimate. However, these two mass estimates tend to disagree 
with the caustic results. We measured the presence of substructure in all clusters of 
the sample and found that clusters with substructure have virial masses higher than 
those derived from Tx ■ This trend is not seen when comparing the caustic and X-ray 
masses. That happens because the caustic mass is estimated directly from the mass 
profile, so it is less affected by substructure. 

Key words: surveys - galaxies: clusters: general - galaxies: kinematics and dynamics. 



1 INTRODUCTION 

Galaxy clusters represent the largest and latest systems to 
form under the influence of their own gravity. The study of 
the evolution of the large scale structure through the clus- 
ter mass function and how that varies with time, is one of 
the most important tools in unveiling the global properties 
of the Universe. The way clusters evolve strongly depend 
on several cosmological parameters like Qm, ^x, and as 



dEke et al.l[l9 98: BahcaU ot al. 1997 , |2003|; iDonahue erall 
1 19981 : iReiciiart et al.. .199a : .Blanchard and Bartlet^ Il998f) . 
The main difficulty encountered in these studies, lies in the 
exact estimate of the cluster mass. A variety of techniques 
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are in use today, but accurate estimates for large data sets 
are still impractical, as all methods are expensive from the 
observational viewpoint. 

A simple way to estimate mass for a large sample of 
clusters is to find an unbiased photometric proxy such as 
X-ray luminosity (La')) richness (Ngai), or optical luminos- 
ity (Lopt). Lx and Lopt correlate well with mass, though 
with larger scatter than relations involving X-ray temper- 
ature (Tx) and velocity dispersions (up). The measure- 
ment and calibration of these relations are essential for re- 
liably establ i shing the cluster mass function. As pointed 
out by IVoitI (|2005h . underestimation of the scatter in the 
mass— observable relation leads to overestimation of erg. 
Furthermore, contradictory results appear when one re- 
lates mass to cluster properties measured at different wave- 
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lengths. This maybe reflecting the selection function of the 
cluster catalogs created from data in different wavelength 
regimes, or the physical processes affecting the cluster galaxy 
population and the intra-cluster gas. 

The study of cosmology using galaxy clusters has as a 
key ingredient the proper recognition of the biases present 
in multi-wavelength surveys. Some of these issues have 
been addressed in the recent literature, either by compar- 
ing X-ray and optical cluster catalogs or c onducting joint 
X-ray/opti c al surveys of galaxy clusters (iDonahue et al.l 
2OO1I. I2OO2I: iGilbank ct al.' 2004': 'Popesso et al."2004l. 12001 



Lopes et al.) [2006: Dai ct al. 2007; Gal ct al. 2008). Some 



of these works investigate the complementary properties of 
clus ters selected in other wavelength regimes (optical or X- 
rav:lponahue et a.l.ll2002|:lGilbank et al.|l20oi:|Popesso et all 



I2OO4I : iLopes et al.ll2006l : lDai et ahlbOOTl - lGal et al.ll2008l ) 



Information on the distribution of gala^xies, hot gas 
and dark matter within clusters can play an important 
role to understand which observational biases contribute to 
increasing the scatter of the scaling relations. The pres- 
ence of substructure is a clear sign of incomplete re- 
laxation in a cluster and can be seen from the X-ray 
emission from the intra-cluster medium (jjones fc FormanI 



ies (IPinknev et al.l 


1996|: 


Mathiesen. Evrard and Mohdl 19991: 


iKolokotronis et al 


2001 


). Recent works have investigated 



the possibility that substructures in clusters increase the 
scatter of t he scaling relations, but the results have been con- 
tradictory. ISmith et al.l (|2005l ) find that the observed scal- 
ing relations for relaxed an d not relaxed clusters are sig- 
nificantly different. However. lO'Hara et al.l (|2006l ) reach the 
opposite concl usion, basing t heir a nalysis on simulated and 
observed data. iLopes et al.l l|2006l ) find that the Tx-Ng^is 
(or Tx-Lopt) relation is severely affected by the presence of 
substructures, although the relations involving Lx are not. 

The main goals of this paper are: (i) collect higher 
quality (SDSS) data for the NoSOCS supplemental cata- 
log to redshift 2: ~ 0.5 (ii) define a low-redshift subset of 
the NoSOCS sample (z ^ 0.1); (iii) derive unbiased mea- 
surements of velocity dispersion and masses for the clusters 
in this low— z sample; (iv) compare optical and X-ray mass 
estimates, allowing checks for possible biases in the mass 
estimates. This paper gives a description of the estimation 
of cluster properties for the NoSOCS catalog using SDSS. 
New photometric redshifts, richness and optical luminosity 
are derived for the full sample. Velocity dispersion, mass, 
and physical radius are obtained for the low-z clusters. This 
data set will be used for future cluster studies. 

This work relies on the su pplement version of NoSOCS 
l|Lopesll2003 : lLopes et al]l2004l '). which contains 9,956 cluster 
candidates over 2,700 square degrees. This catalog, based 
on the Digitized Second Palomar Observatory Sky Survey 
(DPOSS, Djorgovski ct al. 2003), includes clusters up to z ~ 
0.50 and is the largest catalog to this redshift limit to date. 
Here, we discus s the use of highe r quality photometric data 
from the SDSS (|York et al.ll2000D to assess NoSOCS cluster 
properties. We also used spectroscopic data from the SDSS 
to study the dynamics of the low redshift systems (z ^ 0.10). 

The presentation of the paper is as follows: In §2 we 
briefly describe the NoSOCS. Selection of SDSS data for the 
NoSOCS clusters is described in §3, where we also discuss 
the estimates of richness and optical luminosity, as well as 



the effect of different k— corrections and evolution in M* to 
richness. In § 4 we detail the way we characterize the proper- 
ties of the low-z systems (velocity dispersion, characteristic 
radius, and mass). Speciflc issues regarding the method for 
interloper rejection, as well as the radial cutoff and magni- 
tude limit of the spectroscopic survey are also discussed. We 
estimate X-ray luminosity for a subset of these local clusters 
and discuss the estimation of mass for a subset of systems 
with accurate Tx available. In §, 5 we present the substruc- 
ture tests employed for the clusters at low-z. We compare 
optical and X-ray mass estimates in § 6 (as well as opti- 
cal masses from the caustic method and the virial analysis) . 
We summarize our main conclusions in § 7. Throughout this 
work we assumed a cosmology with Q.m =0.3, SI a = 0.7 and 
Ho = 100 h km s'^ Mpc"\ with h set to 0.7. 



2 THE NORTHERN SKY OPTICAL CLUSTER 
SURVEY (NOSOCS) AND CIRS 

The work presented here is based on the supplement ver- 
sion of the N orthern Sky Optic al Cluster S urvey (NoSOCS , 
Lopes 200i iLopes et al.ll2004l ). NoSOCS (|Gal et al.|[200ol . 
2003i . i2008i ') is a catalog of galaxy clusters constructed from 
the digitized version of the Second Palomar Obser vatory Sky 
Survey (POSS-II; DPOSS. loiorgovski et allbOOSl ). This cat- 
alog is derived from high-latitude flelds |fe| > 30°, covering 
~ 11,000 deg^ and containing ~ 15,500 cluster candidates. 
Its construction is limited to r = 19.5, where the star/galaxy 
separation is reliable and the photometric errors are accu- 
rate enough to use the g — r color as a redshift indicator. Ob- 
je ct classification and ph oto metric calibration are described 
in lOdewahn et all Ho^ and lCal et all (120041'). respective ly. 
The supplement version of this survey ( Lopes et al.ll200j ) is 
a deeper cluster catalog (z 0.5) within a smaller area of 
DPOSS (~ 2,700 square degrees), containing 9,956 cluster 
candidates. It covers a smaller area than NoSOCS because 
the search is restricted to the best DPOSS plates, selected 
according to seeing and limiting magnitude {r = 21.0). 
The selection is further restricted to high galactic latitude 
{\h\ > 50°), where stellar contamination is modest and 
nearly uniform. As this catalog is based on deeper galaxy 
catalogs, the larger photometric errors prevent the use of 
colors (from the plate material) for redshift estimates, which 
instead are simply based on a magnitude-redshift relation. 

In section 4 we derive cluster properties for low redshift 
(z ^ 0.1) clusters. In this redshift range NoSOCS comprises 
only poor systems (with a < 700 km/s). That results from 
the small volume probed by the NoSOCS supplement within 
the SDSS spectroscopic survey (at z ^ 0.1). Hence, we de- 
cided to complement o ur sample with the one in cluding mas- 
sive systems studied bv lRines fc Diaferiol(|2006h . Their sam- 
ple comprises 74 X-ray selected clusters at z ^ 0.1 within 
the SDSS DR4 (see their Table 1). This set is named Cluster 
Infall Regions in SDSS (CIRS). That is an optimal data set 
to extend our sample, as it contains more massive clusters 
than this work, being also based on SDSS data and limited 
to the same redshift range. 

Seven clusters of their sample are common to ours; three 
are excluded for being too nearby (z < 0.01) and covering a 
very large angular area. SDSS data for the remaining sixty- 
four CIRS clusters were extracted in the same way as we did 
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for the NoSOCS (section 3). The only difference is the use 
of the SDSS DR6, instead of DR5 (CIRS was incorporated 
into this work when DR6 was available). We refined the 
spectroscopic redshift estimates of these clusters and elimi- 
nated seven systems not well sampled within 2.5 h~^ Mpc 
(inspection performed by PAAL). Another cluster (Zwl665) 
was eliminated for having a large difference between the new 
and the CIRS redshift {\znEw — zcirs\ = 0.0163) and be- 
cause anejD = 0.0139 is too low to be included. The final 
CIRS sample used here contains fifty-six clusters. 



3 SELECTION OF SLOAN DIGITAL SKY 

SURVEY (SDSS) DATA FOR THE NOSOCS 
CLUSTERS 

The photometric and spectroscopic data used in this paper, 
were taken from the fifth release of the SDSS IjYork et al.l 
I2OOOI ). We have selected only objects from the "Galaxy" view 
(so that only PRIMARY objects are allowed) to avoid duplicate 
observations. Standard flags for clean photometry are also 
used. As we select imaging and spectra, a joined query of the 
Galaxy and SpecObj views is performed. All the magnitudes 
retrieved from SDSS are de-reddened model magnitudes. 

Before selecting photometric and spectroscopic SDSS 
data in the region of each NoSOCS cluster, we decided to 
cle an the catal o g in a more rigorous way than was done 
bv lLopes et al.l 1 2004) where we eliminated "double" clus- 
ters by simply excluding objects within 60" of each other 
and located in the overlapping regions of different DPOSS 
plates. These " double" clusters a re candidates with small 
separation fsee iLopes et al.| |2004 ). The originally cleaned 
catalog contains 9,956 cluster candidates. A second clean- 
ing of the NoSOCS catalog was done here, eliminating all 
dusters within 0.75 h"^ Mpc (1.07 Mpc for h = 0.7) of each 
other and a photometric redshift difference of 5C 0.05 (the 
photometric redshift accuracy of the supplemental NoSOCS 
catalog). When two clusters match this criterion, we keep the 
richer cluster. From the 9,956 original systems 358 (3.6%) 
objects are eliminated, leaving 9,598 clusters. 

Considering the higher quality of the SDSS data com- 
pared to DPOSS, and the availability of many colors, we de- 
cided to obtain new photometric redshift estimates (zphoto) 
for all the NoSOCS systems found within SDSS. Thus, ini- 
tially, for each NoSOCS cluster we selected SDSS data in cir- 
cular areas with a radius of 3.0 Mpc -I- 0.5 degrees around the 
cluster center, for which we still consider the original photo- 
metric redshifts from the plate data. We found 7,822 of the 
NoSOCS clusters within the SDSS area. The limiting magni- 
tude of the selected data is r = 21.0, the star/galaxy separa- 
tion limit in SDSS. We have also retrieved data for fifty ran- 
domly selected blank fields of one square degree each. The 
data extraction for these fields is done in the same way as 
for the cluster areas. These control fields are used for back- 
ground subtraction performed when estimating photometric 
redshifts, richness (Ngais) and optical luminosity (Lopt)- 

To obtain new photometric redshifts we adopted an em- 
pirical approach based on the measurement of the mean r 
band magnitude and median (g-i) and (r-i) colors. T he pro- 
cedure for estimating Zphoto is described in section 4 of lLopesI 
l|2007l ). 

With the more accurate photometric redshift estimates 



for each cluster, we repeated the selection of the galaxy data. 
This time we considered regions of 10.0 Mpc x 10.0 Mpc. 
Then, an inspection of finding charts for all cluster areas 
was performed by PAAL. We wanted to exclude clusters 
that were not well sampled, either because they were too 
close to the edges of the SDSS or had a bright (masked) 
object near the center of the field. After this inspection we 
kept 7,550 clusters (96.5%). Finally, we obtain a refined po- 
sition for each cluster. This process of re-centering considers 
the luminosity- weighted coordinat es of all galaxies within 
0.30 h"^ Mpc of each cluster fsee iLopes et al.1 120061 ). This 
aperture is chosen to be large enough to improve the cen- 
troid estimate avoiding projection effects that would affect 
the new coordinates. Photometric redshift estimates from 
the new centroids are also determined for all objects, as the 
Zphoto may be sensitive to this parameter. Considering the 
more accurate redshifts, new Ngais and Lopt estimates are 
obtained for each cluster, as will be discussed in the next 
section. 

As Zphoto is greatly improved from the original value in 
NoSOCS we performed a final cleaning of the cluster cata- 
log, to eliminate "double" clusters. The procedure is analo- 
gous to the one described in the beginning of this section. 
Based on the original coordinates and the associated Zphoto, 
we eliminated clusters that are still within 0.75 h~^ Mpc of 
each other and have a redshift offset of Az ^ 0.025 (tak- 
ing in account the improved accuracy in Zphoto)- Again, we 
keep the richer cluster of two coincident ones. One hundred 
thirty-six (1.8%) of the 7,550 are eliminated. The distribu- 
tion of photometric redshifts for the 7,414 clusters is shown 
in Figure 1, while the richness distribution is exhibited in 
Figure 2. In agreement to the selection function of the cata- 
log, we see that at high redshifts the sample is preferentially 
composed of rich clusters. 

3.1 Richness and optical luminosity 

The procedure to estimate the richness of a cluster depends 
on the cluster centroid, radius adopte d, and the redshift . 
A detailed description can be found in lLopes et al. I (|2006l ). 
Here, we highlight the main points. Our richness measure, 
Ngais, is the number of galaxies at m* — 1 ^5 rrir ^ m* + 
2 within a given aperture, where m* is the characteristic 
apparent magnitude of the cluster luminosity function. For 
this work we consider the bright end values o f the Schechter 
cluste r luminosity function (LF) obtained bv lPopesso et all 
(200(1) , where their fits are optimized for bright and faint 
magnitudes simultaneously. For the bright end, they found 
Q = -1.09 and M* = -20.94 within _R200. These values of 
M* and a are consistent with previous results. The value 
of M* is converted to the same cosmology used here and to 
the proper value at 2: = (taking the mean redshift of their 
sample as z — 0.1). Then, for different redshifts we adopted 
a mild evolutionary correction to the val ue of M*, given by 
M*(z ) = Af*(0) -I- Qz, with Q = -1.4 (|Yee fc Lopez-Cruj 
Il999h . 

As we do not have an estimate of a radius that scales 
with mass (such as R200) we adopt the optimal aperture 
of 0.50 h"^ Mpc (0.71 Mpc for h = 0.7). This choice re- 
sults in the smallest errors in the richness estimates and 
in the most accurate scaling relations (between optical and 
X-ray properties). More details are found in iLopes et al.l 
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Figure 1. Photometric redshift distributions of the 7,414 NoSOCS clusters for which we extract data from SDSS (top). The four bottom 
panels show the z^i^oto distributions in different richness bins (indicated on each panel). 



l|2006l) (similar arguments are given bv lPopesso et al.ll2004 
using the RASS-SDSS sample), k— corrections are ap- 
plied to the magnitude of each galaxy, as explained in 
items (i) and (iii) below. These corrections are obtained 
thro ugh the convolution of spec t ral e nergy distributions 
from IColeman. Wu. fc Weedriianl (|l980l ) with the SDSS r 
filter. The steps to estimate richness are: 

(i) We use the photometric redshift of each cluster to 
convert an absolute magnitude M* to m* and to calculate 



the apparent radius (in seconds of arc) for a fixed aperture 
of 0.50 h^^ Mpc. The photometric redshift is also used to 
compute the k— correction typical of elliptical and late-type 
galaxies (Sbc) . These are named "ke" and "ks" , respectively. 
Since we count only galaxies with m* — 1 ^ m,. m* + 2, 
initially we select all of them within 0.50 h~^ Mpc of the 
cluster center at m* — 1 + fcs ^ rrir ^ m* + 2 + ke. The 
k— corrections are applied to individual galaxies at a later 
stage, so these limits guarantee that we select all galaxies 
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Figure 2. Richness distributions of the 7,111 NoSOCS elusters for which we extract data from SDSS (top). The four bottom panels 
show the richness distributions in different Zphoto bins (indicated on each panel). 



that can fall within — 1 ^ rrir ^ m% + 2. Nciu is the 
number of galaxies selected in the cluster region. 

(ii) We use the fifty blank fields of 1 square degree to esti- 
mate the background counts to be subtracted from the clus- 
ter counts. Galaxies are selected within the same magnitude 
range as used for computing Nciu, and the counts arc scaled 
to the same area of the cluster. The median counts from all 
fifty boxes give the background estimate, ^bkg\ (we actually 
exclude boxes with counts outside the range ±3cr of the me- 



dian counts). We adopt the interquartile range (IQR, which 
is the range between the first and third quartiles) as a mea- 
sure of the error in 'i^bkg, which we term a^kg (for normally 
distributed data IQR = 1.35 cr, where a is the standard de- 
viation). The background corrected cluster counts (Nciu - 
Nifeg) is called Ncorr- 

(iii) Next, we apply k— corrections to the galaxy popula- 
tions in each cluster using a bootstrap procedure. In each 
of the one-hundred realizations, we randomly select Ncor-r 
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galaxies from those in the cluster region (Nc!„). We then 
apply a k— correction to the magnitude of each galaxy. An 
elliptical k— correction is applied to X% of the Ncorr galax- 
ies, while a Sbc k— correction is applied to the remaining 
(100-X)%. Then, we use these k— corrected magnitudes to 
count the number of galaxies at m* — 1 ^ rrir ^ m* + 2. 
The final richness estimate Ngais is given by the median 
counts from the one-hundred realizations. We considered a 
variable fraction of galaxy type at different redshift ranges. 
At 2; ^ 0.f5 we assumed clusters are composed of 80% el- 
lipticals; 50% at 0.15 < z ^ 0.30 and 30% aX z > 0.30. 
Different values for these fractions have little effect on the 
final richness estimates. The richness error from the boot- 
strap procedure alone is given by (Jboot = IQR and the total 
error is the combination of this error and the background 
contribution, so that crjVg^,j, = y^CTb^J^Jt^'-i-'aftfe^. From now 
on, we simply call the richness error Ngau-err- 

(iv) If the cluster is too nearby or too distant, either the 
bright (m* — 1-1- ks) or the faint (m* -f 2 -f- ke) magnitude 
limit will exceed one of the survey limits (10.0 ^ 21.0). 

Here, we apply one of the following correction factors to the 
richness estimate; 



71 



72 



J m^, — 1 ^ ' 

Iw'^^ $(m)dm 

Jm^ — 1 ^ ' 
Jm:z — 1 ■ ' 



(1) 



(2) 



We call 7i and 72 the low and high magnitude limit correc- 
tion factors. Whenever necessary, one of the above factors is 
multiplied by Ngais and Ngais-err- 

The photometric redshift and richness distributions are 
illustrated in Figures 1 and 2. Note that the distributions 
use richness estimates without evolution in the value of M*. 
Such evolution would make the richness estimate smaller and 
would cause the distributions to be shifted to the left. These 
differences affect only clusters at higher redshifts (z > 0.25). 
Although the listed values in the tables consider evolution 
in M*, the mass estimates and scaling relations, obtained 
only at 2: ^ 0.10, are not affected by this correction. 

The DR5 of SDSS also provides photometric redshifts, 
galaxy types and k— corrections for all galaxies obtained 
through a template fitting procedure. Thus, we can test if 
the richness estimates are sensitive to the k— correction used, 
if individual or based on a statistical approach. We consider 
the individual k— correction and zphot of each galaxy to 
obtain new estimates of richness. The procedure is analo- 
gous to that used above, except for the k— correction that 
is applied to each galaxy (see item iii above). In all the 
one-hundred realizations, when selecting Yicorr galaxies we 
apply the k— correction given by SDSS to each galaxy, avoid- 
ing the need to statistically consider different galaxy types. 
However, as we can see from Figure 3, the richness estimates 
obtained in both cases are indeed very similar. 

To estimate the optical luminosity of each cluster, we 
proceed as follows. We generate the magnitude distribution 
(assuming a bin of Amag = 0.2) of the Ncorr galaxies found 
within the cluster regions, at each magnitude bin. Then, the 
corrected optical luminosity is given by 



L = ^ivao-°- 



4r^ 



(3) 



where Ni represents the corrected counts at each magni- 
tude bin Ti. The sum is performed over all n magnitude 
bins. This luminosity is converted to an apparent magni- 
tude associated to the cluster, from which we subtract an 
elliptical k— correction, rridu = —2.5log{L) — ke. We then 
consider the distance modulus (DM) for the cluster red- 
shift and transform this value to an absolute magnitude, 
Mciu = rriciu — DM. Finally, considering the absolute mag- 
nitude of the Sun in t he r band of SDSS (M^ = 4.64; 
iBlanton fc Roweid [20071 ) ■ we obtain the optical luminosity 
in solar units as L^pt = 10"°"'*^='""*^o' . The error in the 
luminosity is given by 



AL : 



(4) 



where \/77i and aibkg are the errors in the cluster corrected 
counts and background in the i-th magnitude bin. The term 
in brackets represents the combined error in each bin. If 
the survey limits are exceeded, one of the correction factors 
below is multiplied by L (equation 3) and AL (equation 4). 



Ti 



(5) 



(6) 



4 PROPERTIES OF LOW REDSHIFT 
CLUSTERS 

In this section, we describe the procedures adopted for es- 
timating cluster properties such as velocity dispersion (a), 
characteristic radius (-R500 and ii20o) and mass (M500 and 
Af20o)- These parameters are measured only for clusters at 
z ^ 0.10. We use the SDSS spectroscopic survey to derive 
these estimates. This survey is complete to z ~ 0.10, where 
the peak of the redshift distribution is located. At higher 
redshifts we start missing galaxies fainter than Af * -I- 1, and 
that may bias the dynamical analysis. Details regarding this 
choice are given below. It is also important to stress that 
some of the more recent cluster catalogs from SDSS are not 
always well sampled at z ^ 0.10. For instance, the MaxBCG 
catalog was constructed at 0.10 ^ z ^ 0.30. 

The first step to determine velocity dispersion, char- 
acteristic radius, and mass is the identification of the clus- 
ter spectroscopic redshift and velocity limits, as well as the 
selection of cluster members from the spectroscopic sam- 
ple. We start by identifying groups in redshift space and 
thus obtain a spectroscopic reds hift for each cluster . We 
use the gap-technique described in lKatgert et al.l (|l996l ) and 
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Figure 3. Comparison between richness estimates obtained with 
a statistical k— correction (horizontal axis) and using the values 
obtained for each galaxy in SDSS (vertical axis). The solid line 
(red in the electronic edition) indicates the Y = X result. 

lOlsen et all (|2005[ ) to identify gaps in the redshift distribu- 
tion. Instead of adopting a fixed gap, such as 1000 km 
or 0. 005(1 + z), we cons idered a variable gap, called density 
gap (jAdami et al.lll99"^ . The density gap size is given by 
the expression Az — 500(1 + exp( — (A'^ — 6)/33))/c, where 
N is the number of galaxies found in the redshift survey of 
a cluster (|Adami et al.|[l998l ). and c is the speed of light in 
km/ s. Here, we select all galaxies within 0.50 h~^ Mpc of the 
cluster center. When performing galaxy selection, the only 
requirement is that the galaxy is brighter than r = 18.0 (ap- 
proximately the magnitude limit of the SDSS spectroscopic 
survey). In the region of each cluster, groups in z-space are 
identified through the application of the gap-technique. The 
redshift of each g roup is given by the biweight estimate 
l|Beers et al.l |l990!). 

The next step is to assess the significance of each of 
these groups. For that we consider the full spectroscopic 
survey of SDSS (considering the same limit at r = 18.0). 
For a given cluster, for each group identified within it, we 
draw 1000 sets of galaxies. These sets have the same number 
of galaxies as in the cluster region where the group was 
identified. The gap-technique is applied exactly as before, 
and then we check the probability of finding groups with 
at least the same number of galaxies at the redshift of the 
original group. We consider only groups that are significant 
at the 95% level. 

From all the significant groups (sometimes there is only 
one), we select the one that has the smallest redshift dif- 
ference to the photometric value of the cluster. We also 
require this difference {\zphoto — Zspecj) to be smaller than 
0.03(1 + Zspec), where Zp^oto is the photometric redshift esti- 
mate of the cluster (§ 3.1) and Zspec the spectroscopic value 
of the group. Finally we also require that the group have at 



least three member galaxies within the 0.50 h~^ Mpc aper- 
ture (most have many more). During this procedure, we also 
obtain the velocity (or redshift) limits [vio and vm) of the 
cluster within this radius. This procedure is applied to all 
754 clusters (out of the original 7,414) with Zphoto ^ 0.133. 
From those we find 179 clusters with z^pec ^ 0.10 (hav- 
ing at least 3 galaxies within 0.50 h"'^ Mpc). From now on, 
we will refer to z^pec as Zduster- We noticed that for two 
systems the velocity dispersion estimates (see below) were 
artificially high because of the influence of a projected clus- 
ter along the line-of-sight. We found that a refined redshift 
estimate (within 0.30 h"^ Mpc, instead of 0.50 h"^ Mpc) 
helps solve this problem. 

Next we apply an algorithm for interloper rejection to 
arrive at a final list of cluster members, which will be used 
to estimate the cluster properties. The procedure we adopt 
is similar to that of [Fadda ct al. (1996) and is called the 
"shifting gapper" technique. It works through the applica- 
tion of the gap-technique in radial bins from the cluster 
center. The bin size is 0.42 h"^ Mpc (0.60 Mpc for h = 
0.7) or larger t o force th e selection of at least 15 galaxies 
(consistent with lFadda et al. 1996). Galaxies not associated 
with the main body of the cluster are eliminated. This pro- 
cedure is repeated until the number of cluster members is 
stable (no more galaxies are rejected as interlopers). This 
method makes no hypotheses about the dynamical status of 
the cluster, while other procedures are based on physical as- 
sump tions about the cluster mass profile (see lWoitak et al] 
I2OO7I ). The way we apply this technique has many details, 
which are described below. 

(i) We select all galaxies within 2.5 h"^ Mpc (3.57 Mpc 
for h = 0.7) from the cluster center and with \ cz—czciuster\ ^ 
4000 km . The aperture of 2.5 h~^ Mpc is typically larger 
than the virial radius of clusters, but is still small enough 
to avoid the influence of large scale structure on the cluster 
analysis. More details about this choice are given below. 

(ii) Before starting the shifting gapper procedure, we take 
the maximum velocity offset and impose the velocity lim- 
its to be symmetric. The velocity offsets within 0.50 h~^ 
Mpc are \vio ~ Vdusterl and \vhi — Vduster], where v duster = 
czciuster- The maximum offset between these two is named 
Vdiff and the new symmetric velocity limits are defined 

as Vlo = Vduster — Vdiff aud Vhi = Vcluster + Vdiff- If 

Vdiff < 300(1 -I- Zduster) km we set it to this value. 
We have also noticed that the factor "500" in the density 
gap expression (Az — 500(1 + exp(— (A*' — 6)/33))/c) may 
not be suitable when running the shifting gapper for low and 
high mass systems in radial bins. Thus for each radial bin 
we define the parameter facgap = \vio — Vhi\/W and modify 
the expression to be Az = facgap{l + exp{—{N — 6)/33))/c. 
This parameter facgap assumes a minimum value of 300. 

(iii) When we select galaxies within a radial bin of 0.42 
h~^ Mpc or larger (in the case less than fifteen galaxies are 
found), if the radial offset between two consecutive galaxies 
is greater than 0.7 h~^ Mpc, the procedure is halted. Then 
all galaxies with a radial distance from the cluster center 
greater or equal to this galaxy's radial distance are elimi- 
nated as interlopers. 

(iv) When starting the procedure, in the first radial bin, 
we define vio, v^i and facgap as above. Then we run 
the gap technique for the galaxies of that radial bin. If 
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at least one group is identified in redshift space, we sort 
these in velocity offset from the cluster. We consider the 
first group (smallest velocity offset) valid if the offset is 
^ c(0. 001(1 + Zciuster))- We thcu consider all the galaxies 
from this group as cluster members, except those outside 
the velocity limits, vio and vm- If the first group offset is 
larger than c(0. 001(1 + Zduster)), or if there are no groups 
found in this radial bin, we still keep all galaxies within vio 
and Vhi- If the first group is considered as valid, we also take 
the velocity limits of this group to update the values of vio 
and Vhi, exactly as above, keeping these values symmetric 
and varying with radial distance. For the next radial bin, 
these new values of vio and Vhi will be considered and so on. 

After the interloper removal we kept one-hundred and 
twenty-seven clusters (out of one-hundred and seventy-nine) 
with at least ten member galaxies selected. Figure 4 illus- 
trates the result of the cluster member selection for these 127 
NoSOCS systems. On each panel one galaxy cluster is rep- 
resented, with filled squares showing cluster members and 
open circles galaxies rejected as interlopers. 

We then perform a virial analysis of all the one-hundred 
and twenty-seven clusters with at least ten galaxies (trying 
to use systems with less than ten galaxies leads to increased 
scatter in the scaling relations). Th e proc e dure is analo- 
TOUstowhat is d one bylGirardi et al ]Jl998l ). |Po"pesso et al.l 



ROUS to wna t is d one bvHjirardi et ai.i li^^ai i. iJr'opesso et ai.l 
l|2005l . bOOTl ) and iBiviano et all l|20o'g 'l. First, we compute 
the line-of-sight velocity dispersion (crp) of all cluster mem- 
bers within the aperture Ra, which is the radial offset of 
the most distant cluster member. Ra is usually close to 2.5 
h~^ Mpc (the maximum radius adopted for this work), but 
could be much smaller if the outermost galaxies are con- 
sidered interlopers. The robust velocity dispersion estimate 
(crp) is given by the gapper or biweight estimator, depending 
on whether the number of objects availab le is < 15 (gap- 
per) or ^ 15 (biweight; Beers et al.l [l990l '). Following the 
prescriptions of lDanese et al.1 (I198G| )7 the velocity dispersion 
is corrected for velocity errors. We als o obtain a n estimate 
of the projected "virial radius" Rpv fjGirardi et al.lll99S) . 
A first estimate of the virial mass is given by (equation 5 of 
iGirardi et al]ll998l ) 



Mv 



2G ' 



(7) 



where G is the gravitational constant and 37r/2 is the de- 
projection factor. 

Next we apply the su rface pressure term correc- 
tion to the mass estimate l)The fc White! [l986'l. To do 
that we first estimate the Navarro et al. (1997, NFW 
from now on) concentration paramete r, using the relation 
c = 4 X {M/ Mkbm)~°-^°'^ (see also iPopesso et al.1 [2OO7I: 



iBiviano glalj [looi) . The slope is taken from iDolag et al.l 



(2004) and the nor r naliza tion [Mkbm « 2 x 10^^ Mq) is 
from iKatgert et al.l (|2004l 'l. To apply the surface pressure 
term we also need an estimate of the R200 radius, which (as 
a first guess) is taken from the definition of Carlberg et al. 
(1997; equation 8 o f that paper). The c orrection (computed 
as in equation 14 of lGirardi et al.l|l998l ) is applied consider- 
ing isotropic orbits. In particular in our analysis, the term 
associated to the presence of v elocity anisotropies is set to 
0.3, as found by iGirardi et al.l (Il998{) at R„i,.. We noticed 
that for different velocity dispersion profiles, and at other 



radii, different values should be adopted for this term. That 
may lead to differences of up to ~ 25% in the value of the 
surface pressure correction, but the effect on the mass is less 
than 5%. Since the effect is much smaller than the typical 
error in mass, and the study of velocity anisotropies is be- 
yond the scope of this paper; we consider the value of 0.3 
suitable. Note also that the most important issue here is to 
apply the surface pressure correction, which is sometimes 
neglected in the literature. 

After applying this correction, we obtain a refined 
estimate of R200 considering the virial mass density. If 
Mv is the virial mass (after the surface pressure correc- 
tion) in a volume of radius Ra, R200 is then defined as 
R200 ~ -Ra[pv/(200/9c(2))]^^^'*. In this expression pv = 
3Mv /{4:TvR\) and Pc{z) is the critical density at redshift 
z. The exponent in thi s equation is the on e describing an 
NFW profile near R200 (iKatgert et al.ll2004l ). If we use five- 
hundred instead of two-hundred on this equation, we obtain 
an estimate of Rsoo- Next, assuming an NFW profile we ob- 
tain M200 (or M500) from the interpolation (most cases) or 
extrapolation of the virial mass Mv from Ra to R200 (or 
Rsoo). Then, we use the definition of M200 (or M500) and a 
final estimate of R200 (or R500) is derived. Thi s proc edure 
is analogous to wha t is done by iPopesso et al.l (|2007| ) and 
IBiviano et al.l (I2OO6I ). We assume the percentage errors for 
the mass estimates are th e same before and a fter the surface 
pressure correction (as in lGirardi et al.lll998l ). 

Table 1 lists the one-hundred and twenty-seven 
NoSOCS clusters and their main characteristics. The cluster 
name is in column 1; coordinates are shown in columns 2 and 
3; redshift in 4; richness and its error in column 5; optical 
luminosity and its error in column 6; and the value of the 
j3 and A parameters obtained from the substructure tests 
and their significance (see section 5) are given in columns 
7-10. Table 2 lists the velocity dispersion, number of cluster 
members used for that estimate, number of cluster members 
within R200, and the characteristic radii and masses of these 
clusters, R500, M500, R200, and M2oo- 

Table 3 lists the main characteristics of the CIRS sys- 
tems. The same parameters, as given in Table 1, are shown 
here plus the original redshift listed in the CIRS (fourth col- 
umn). Table 4 is analogous to Table 2, but for the fifty-six 
CIRS clusters. Phase-space diagrams for these fifty-six clus- 
ters are shown in Figure 5. Note that the one-hundred and 
twenty-seven NoSOCS clusters have velocity dispersion esti- 
mates of 100 < (T < 700 km/s, while the CIRS systems have 
200 < a < 900 km/s (with only 23% of objects with a < 400 
km/s). We noticed that three CIRS systems still have a large 
difference between the original and new redshifts, but not as 
large as Zwl665. For one of these three clusters, we find that 
a new redshift estimate within 0.3 h~^ Mpc minimizes the 
influence of a neighbor system (as done for two NoSOCS 
clusters; see above). These three clusters are outliers in the 
scaling relations (paper II). However, they are still listed in 
Tables 3 and 4. These three clusters are Abell 1035B, Abell 
1291 A, and Abell 1291B (see Figure 5). 

In the next three subsections, we discuss a few issues 
that could affect the estimates of cluster properties. We fo- 
cus on the interloper rejection and radius and luminosity 
limits of the spectroscopic survey. 
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Figure 4. Phase-space diagrams of 15 NoSOCS clusters shown as examples. The velocity and radial offsets are with respect to the 
cluster center. We apply a shifting gapper procedure for the selection of cluster members (filled squares) and exclusion of interlopers 
(open circles). Similar diagrams for all NoSOCS clusters are available in the electronic edition of the MNRAS. 



4.1 Testing the rejection of interlopers 

In the literature, we find several different approaches for in- 
terloper removal in the dynamical modeling of galaxy clus- 
ters. For a recent comparison of the performance of several 
different methods a pplied to A'^— body c osmological simula- 
tions see the work of lWoitak et al.l ()2007l ) . As noted by these 
authors, differences in mass estimates may be explained by 
the number of interlopers a given method selects. That could 



also be an explanation for discrepant estimates to other 
methods based on X-ray observations or lensing analysis. 
Here, we do not intend to perform a detailed comparison. 
We use a method (shifting gapper) that has the following 
analysis advantages: use of the combined information of po- 
sition and velocity; and is independent of any hypotheses 
regarding the dynamical status of the cluster. However, it is 
important to mention that even slight modifications in the 
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Table 1. Main properties of the 127 NoSOCS clusters at z sC 0.01 and with at least 10 member galaxies. 



name 


ra 
(J2000) 


dec 
(J2000) 


Z0 


N I 

gals 


Lopt 
(lO^^Lo) 


P 




A 




NSCS J121847+484410 
NSCS J011502+002441 
NSCS J100242+324218 


184.67720 
18.75351 
150.67372 


48.72197 
0.37955 
32.70033 


0.0448 
0.0449 
0.0505 


10.0 ± 3.0 
66.0 ± 14.5 
29.0 ± 7.8 


0.229 ± 0.093 
1.256 ± 0.213 
0.351 ± 0.108 


109.600 
46.200 
-14.400 


0.048 
0.000 
0.460 


6.577 
79.405 
44.065 


0.717 
0.317 
0.571 



Note. - A portion of this table is shown here for guidance regarding its form and content. A full version is available in the 
electronic edition of the MNRAS. 



Table 2. Velocity dispersion, characteristic radii and masses of 127 NoSOCS clusters. 



name ap Nobj-o-p Nobj-R200 R-SOO M500 R200 M200 

(km s-i) (Mpc) (IO^^Mq) (Mpc) (lO^^M©) 



NSCS J121847+484410 179.9 +11'}^ 29 12 0.49 to.og 0.35 to^l 0.67 to gg 0.36 to^t 

NSCS J011502+002441 474.6 tfsl US 78 1.05 to°ol 3.46 ^qaI 1-45 lo'oe 3.63 ttlt 

NSCS J100242+324218 376.1 tfii 78 46 0.78 +°;°^ 1.41 1.07 +°;°g 1.45 Ig'g 

Note. - A portion of this table is shown here for guidance regarding its form and content. A full version is available in 
the electronic edition of the MNRAS. 



procedure described above could lead to different estimates 
of the cluster properties. For instance, had we decided not to 
update the velocity limits and Vhi from one radial bin to 
the next, we would end up with fewer interlopers (see item 
iv of previous section). Then we would keep the values of vio 
and Vhi found within 0.50 h~^ Mpc, instead of making those 
smaller for larger radius. If we had also enlarged these ini- 
tial limits by a 3% factor (to account for uncertainties in the 
choice of these limits) , the number of interlopers would also 
decrease. For some clusters, the velocity dispersion could be 
~ 20% larger and the mass ~ 50%. For comparison, Figure 
6 shows the phase-space diagrams for one cluster as an ex- 
ample. In the top panel, we show the same distribution as in 
Figure 4 and in the bottom we consider the velocity limits 
(within 0.50 h~^ Mpc) to be larger by a factor of 3% and do 
not update those for larger radial bins. As we can see, the 
number of interlopers (open circles) is now smaller. We de- 
cided to keep the approach described in the previous section 
because it is more stable at larger radii (especially for low 
mass systems). The scaling relations (see lLopes et al.ir2008l : 
paper II) obtained with this approach are also found to have 
smaller scatter. 
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4.2 Radial cut-off 

One can consider as physically meaningful only the values of 
velocity dispersion found after considering all cluster galax- 
ies within a sufficiently large radius . This issue has been 
previo usly discussed by 
l|l996l ) and iGirardi et al 



Fadda et all (|l996l ). iGirardi et all 
( 1998l l. Here, we investigate how 



the estimates of ap are affected by the maximum aperture 
(R-max) considered. In Figure 7, we show the comparison of 
ap obtained with Kmax = 2.5 h~ Mpc (horizontal axis) and 
Rmax = 0.5, 1.5 and 3.5 h~^ Mpc (vertical axis). This fig- 
ure clearly shows that a small aperture (0.5 h~^ Mpc) tends 
to overestimate the velocity dispersion. At larger apertures, 
the effect is more subtle. We use Rmax = 2.5 h~^ Mpc be- 
cause it is typically larger than the virial radius of clusters 
and still small enough to avoid the influence of large scale 
structure or nearby clusters. From inspection of the phase- 



Figure 6. Phase-space space diagrams of the cluster NSCS 
J121742-I-034017. In the top panel we show the original results 
(same as in Figure 4), while in the bottom we show the distri- 
bution when considering a less rigorous criteria for interloper re- 
moval (see section 4.1). Although the differences in the final num- 
ber of galaxies rejected (open circles) are small, the larger number 
of cluster members (filled squares) results in larger estimates of 
velocity dispersion and mass. In particular, in a few cases those 
estimates could differ by more than 50%. 



space diagrams of several clusters, we found this aperture 
to be a good choice when working both with low and high 
mass systems. 
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Table 3. Main properties of the 56 CIRS clusters at z sC 0.01 and with at least 10 member galaxies. 



name 


ra 
(J2000) 


dec 
(J2000) 


Z0 

(CIRS) 


Z0 

(New-SDSS) 


N I 


(IOI^Lq) 


P 


0sig 


A 




MKW4 

MKWll 

A0779 


181.10500 
202.38361 
139.92200 


I. 90060 

II. 78920 
33.76300 


0.0204 
0.0228 
0.0233 


0.0201 
0.0229 
0.0231 


27.0 ± 8.1 
37.0 ± 8.2 
33.0 ± 7.9 


0.461 ± 0.168 
0.519 ± 0.121 
0.546 ± 0.120 


-8.400 
-4.700 
6.000 


0.540 
0.402 
0.194 


106.903 
62.810 
107.662 


0.052 
0.140 
0.138 



Note. - A portion of this table is shown here for guidance regarding its form and content. A full version is available in the 
electronic edition of the MNRAS. 



Table 4. Velocity dispersion, characteristic radii and masses of the 56 CIRS clusters. 



name ap J^obj-crp Noi,j-R200 R-SOO M500 R200 M2()0 

(km s-i) (Mpc) (lO^^'M©) (Mpc) (lO^^'M©) 

MKW4 392.4 +fg^o ^35 81 0.88 to.o4 1-97 tofg 1-21 to.oe 2.04 to.to 

MKWll 359.4 tlil 119 55 0.91 +°'°* 2.19 toll 1-25 ^oo4 2.27 

A0779 400.2 135 81 0-89 2.07 +«;^^ 1-23 +°;;^^ "^-^5 ttfv 

Note. - A portion of this table is shown here for guidance regarding its form and content. A full version 
is available in the electronic edition of the MNRAS. 




Figure 7. Comparison of velocity dispersion estimates obtained 
for different maximum apertures to select galaxies from a spec- 
troscopic survey around clusters. In the horizontal axis, we show 
the results considering Timax = 2.5 h~^ Mpc. In the top panel 
the comparison is to the values of ap with Rmaa; = 0.5 h~^ Mpc, 
while with Umax = 1.5 h~^ Mpc and Umax = 3.5 h~^ Mpc are 
shown in the middle and lower panel, respectively. The Y = X 
line is indicated in the three panels. 

4.3 Luminosity limit 

Now we want to see how the magnitude limit of a given spec- 
troscopic survey may aflect the velocity dispersion estimates 
(as well as R200 and M2oo). From inspection of the magni- 
tude and redshift distributions of the SDSS data, we find 
that the survey is approximately complete to r = 17.6 and 
z = 0.1, in agreement to the stated values {rpetro ~ 17.77). 



At z = 0.1 this magnitude limit allows one to sample M* -|- 1 
galaxies. At the bright-end of the luminosity function we 
expect th e SDSS spectros copic sample to be incomplete at 
r < 14.5 (|York et al.ll200ol ). 

To test the effect of missing galaxies at the bright or 
faint end of the luminosity function, we simply cut the SDSS 
data at different magnitude limits. For clarity, we show in 
Figure 8 only the results when considering galaxies brighter 
than r = 16.6. On every panel, the vertical axis has the 
results obtained with the full spectroscopic sample, while 
the horizontal axis has the results for galaxies brighter than 
r=16.6. The upper panel shows the comparison for the ve- 
locity dispersion, while the characteristic radius (R200) and 
mass (M200) are exhibited in the middle and lower panel, 
respectively. When restricting the sample to galaxies at 
r ^ 14.5, the velocity dispersion estimates are not affected 
by missing a few bright objects, even for the low redshift 
clusters. However, when selecting only galaxies at r 16.6 
we notice that many clusters have their estimates affected by 
missing galaxies in the range M* < M < M* 4- 1. The filled 
circles in the figure represent clusters at z 0.065, which 
sample at least M* -I- 2 with the full spectroscopic sample 
(complete to r = 17.6) and M* -|- 1 for the limited sample 
(at r = 16.6). We see that for these clusters the scatter is 
smaller compared to the open circles, indicating that miss- 
ing galaxies at M* -I- 1 < M < M* -\- 2 have only a minor 
effect on the ap estimates (as well as in the values of R200 
and M2oo). 

Using the results shown in Figure 8 we can compute 
the mean and standard deviation of the difference between 
the velocity dispersion estimated with the full spectroscopic 
sample and only galaxies brighter than r — 16.6. We find 
fi — 3.09 km and a = 95.74 km . In terms of mass, 
we have = 0.02 lO^"* M0 and a = 1.55 lO" Mq. Ahhough, 
the offset is consistent to zero we stress the large scatter. If 
we make a cut to consider only clusters with ap ^ 300.00 
km (velocity dispersion of the X axis in Figure 8) we 
have the following: for velocity dispersion fi — 28.89 km 
and a — 85.11 km , while for mass fi — 0.32 lO^** Mq 
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Figure 5. Phase-space diagrams of fifteen CIRS clusters arc shown here. The velocity and radial offsets are respective to the cluster 
center. We apply a shifting gapper procedure for the selection of cluster members (filled squares) and exclusion of interlopers (open 
circles). Similar diagrams for all CIRS clusters are available in the electronic edition of the MNRAS. 



and a = 1.11 10^"* M©. From that wc SCO a small trend for 
underestimation of ap and mass for low mass systems. 

In the ideal case, we would like to have a spectroscopic 
survey as deep as possible. However, when that is not possi- 
ble we can see from Figure 8 that a survey which is complete 
to M* + 1 is sufficient to get unbiased estimates of crp, R200 
and M200- Figure 9 shows the phase-space diagram of the 
most massive NoSOCS cluster studied here, at 2; = 0.0776. 



In the upper panel we see all galaxies from the full spec- 
troscopic survey, while the lower panel has the distribution 
of galaxies with r ^ 16.6. We see that the absence of the 
M* < M < M* + 1 galaxies, in this case, transforms a mas- 
sive cluster into a low richness system. The original estimates 
for this cluster are ap = 693.57 km s~^ and M200 = 9.07 
10" Mq. With the reduced data we find ap = 427.28 km 
and M200 = 3.69 10^* M0. It is thus clear that using 
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Figure 8. Comparison of velocity dispersion, characteristic ra- 
dius, and mass estimates obtained for different cuts in luminos- 
ity. The vertical axis in all panels has the results with the full 
spectroscopic sample, while the horizontal axis has the results 
for galaxies brighter than r = 16.6. The filled circles represent 
clusters at z ^ 0.65, which sample at least M* + 2 with the full 
spectroscopic sample (complete to r = 17.6) and M* -|- 1 for the 
limited sample (at r = 16.6). The Y = X line is indicated in all 
the panels. 
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Figure 9. Velocity diagram of the most massive NoSOCS cluster 
in the sample studied here. In the upper panel, we show the orig- 
inal results obtained with the full spectroscopic survey of SDSS 
(same as in Figure 4). In the bottom panel we show only galaxies 
at r sC 16.6, to see how the interloper rejection, velocity disper- 
sion, and mass estimates are affected. The original values for this 
cluster are ap = 693.57 km s'^ and M200 = 9.07 lO^"* M©. Ex- 
cluding galaxies at r > 16.6, we obtain ap = 427.28 km and 
M200 = 3.69 10^^ Mq. 



galaxy clusters at redshifts where the spectroscopic survey is 
very incomplete will severely affect the final mass estimates. 

4.4 X-ray luminosity from the ROSAT All Sky 
Survey (RASS) 

Following iGal et all l|2008h we used the locations of opti- 
cally selected cluster candidates to measure X-ray fluxes 
and luminosities from RASS. The significance of X-ray 
emission in some of these areas may be too low to iden- 
tify extended sources. However, data from this survey are 
still useful for deriv i ng eit her fluxes or upper limits follow- 
ing |B6hringer|eF^al] (|2000l ). Although the X-ray data from 
RASS are limited, especially for the lower mass systems, 
the knowledge of the cluster positions allows us to esti- 
mate X-ray measurements for a large portion of our sample, 
comprised of 127 NoSOCS clusters plus 56 CIRS systems. 
Note that all CIRS objects had X-ray luminosities previously 
measured (obtained from diflterent X-ray catalogs) . However, 
we decided to re-estimate those in order to have a homoge- 
neous set of estimates. 

The X-ray luminosities Lx are estimated from count 
rates in ROSAT PSPC images. Images and exposure maps 
are retrieved from the ROSAT archive via FTP. We begin 
by estimating the background contribution. That is done 
in three different ways. The first estimate is identical to the 
procedure adopted for the Northern ROSAT All- Sky Galaxy 
Cluster Survey (NORAS. iBohringer et al]|200d ). where the 



background is determined from a ring area centered on the 
source. The inner ring radius is 21' and the outer is 41.3'. 
The annulus area is further divided in 12 sectors, and the 
background is estimated from the median counts of all sec - 
tors. Further details can be found in lBohringer et ah l (|2000l ). 
The second background estimate is based on the selection 
of a hundred boxes of one-hundred square arc minutes over 
the RASS field of 6.5° x 6.5° where the cluster is. In both 
cases, the median and rms (bkgmed and bkgrms) are taken 
(from the 12 sectors or the 100 boxes). The third approach 
is based on the counts obtained in this whole frame of 6.5° x 
6.5°. For this work we estimate the flux for each cluster in 
three different apertures, 0.5 h~^ Mpc, R500 and R2oo- Only 
clusters whose total counts are 3a above the background are 
considered reliable. All others are reported as upper limits. 

To convert the measured total count rate into an un- 
absorbed X-ray flux in the nominal ROSAT energy band 
(0.1-2.4 keV) we use the PIMMS tool available through 
NASA H EASARC. We assume a Raymond- Smith (RS) 
spectrum l|Ravmond fc SmithI 1 19771 ) to represent the hot 
plasma present in the intracluster medium, with a metal- 
licity of 0.4 of the solar value (the results are insensitive 
to small variations in metallicity). The interstellar hydro- 
gen column density along the line-of-sight is taken as the 
weighted average value of iDickev fc LockmanI (|l990l ) (that 
is selected using the 'X-Ray Background Tool", also avail- 
able through NASA HEASARC). The flux in the nominal 
ROSAT energy band is calculated assuming a flxed temper- 
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Figure 10. Comparison of X-ray luminosities obtained for tliree 
different background estimates. In the lower panel the results con- 
sidering the "boxes" background is compared to the "annulus" , in 
the central panel the comparison is between the "frame" and "an- 
nulus" . In the upper panel luminosities obtained with the "frame" 
are compared to the "boxes" . 



ature of 5 keV, which is typical for clusters (|Markevitchl 
1 19981 ). The resulting luminosity is termed Lxs- Then we 
use an iterative proced ure relying on the Ljf-Tx relation 
from iMarkevitchI (|l998l ). From Lxs we find a new tempera- 
ture which is used to recalculate the luminosity based on 
an RS spectrum. The procedure is iterated until conver- 
gence is reached (when the change in temperature is Tx < 
IkeV , comparable to the scatter in the Lx-Tjc relation). 
The procedure typically converges in two or three iterations. 
For both luminosity measures, we apply a k— correction 
l|Bohringer et al.ll2000l ) to derive the X-ray luminosity in the 
rest frame 0.1 - 2.4 keV band. 

Figure 10 shows a comparison of the X-ray luminosi- 
ties obtained for the three different background estimates 
discussed above. As we can see the luminosities obtained 
with the "annulus" tend to be lower than the two other 
cases for the most poor systems. However, for systems with 
LH200 ^ Q _ Q 2 iQ'ii gj.g/g the "annulus" and "boxes" 
are consistent to each other. When the comparison regards 
the "frame" the luminosities with the "annulus" background 
are always smaller. The results with the "frame" are also 
generally higher than those based on the background ob- 
tained with the "boxes". So, we conclude that the "frame" 
background generally provides luminosities that are overes- 
timated compared to the other two backgrounds. The good 
agreement between "annulus" and "boxes" is also seen in 
Figure 11, where we show the ratio of luminosities as a func- 
tion of R200, ^^aia ""f • It is clear that for the richer 
systems {^^aia > 40) the X-ray luminosities are more reli- 
able. 



Figure 11. Ratio of X-ray luminosities obtained with the "boxes" 

and "annulus" backgrounds as a function of R200 (lower panel), 
ii200 

gals 



N^^'l" (middle) and crp (upper panel) 



Table 5 lists the values of Lx (and its associated error) 
for the three backgrounds; "annulus" , "frame" and "boxes" , 
respectively. We also list the X-ray temperature measure 
from BAX (see section 4.5) when available and the interpo- 
lated temperature obtained from the Lx-Tjc relation. The 
last column of this table indicates whether we had a sig- 
nificant detection (SD), an upper limit (UL), or the source 
is too close to the border implying that the measurements 
might be taken with concern (XX). The last two columns are 
obtained with the background estimated from the anrmlus. 



4.5 Cluster masses estimated from Tjc 

We searched BAX [Base de Donnees Amas de Galax- 
ies X, http://bax.ast.obs-mip.fr/) for counterparts of 
the 183 NoSOCS and CIRS clusters used in this work. 
We restricted the search to clusters sd, z < 0.12 with 
X-ray temperature measures available. We found 282 
clusters in BAX, of which 21 are common to our sample. 
Actually, there were twenty-five clusters in common, 
but we rejected four that had no error estimates on 
their temperature. Dat a for these twent y-one systems 
come from Einstein (David e t al ] I1993 I). BeppoSAX 
Mariiii et all [200 4). ASC A (Finoguenov.Arnaud. fc Davidl 

)hashil 



2001 



Ikebe et al.l l2002l: Fukazawa. Makishima. fc Ohash: 
200411 ■ XMM jBelsole et al 



- .20041) and Chandra 

llMarkevitch et all lioocl: lYang et al.l l2004l : IVikhlinin erall 



I2OO6I: iFuiita. Sarazin. fc Sivakofj l2006t l. The temperature 
values for these 21 clusters are listed in the fifth column 
of Table 5. To check the homogeneity of the temperature 
measures we selected from BAX, we show in Figure 12 
the ratio betwe e n th e temperature measures listed in 
iRines fc Diaferid (|2006l ) and BAX. The former is retrieved 
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Table 5. X-ray luminosity and temperature of the 183 NoSOCS plus CIRS clusters. 



name 


Lx (annulus) 


hx (frame) 


Ljf (box) 


Tx (BAX) 


Tx (interp. RASS) 


NOTE 




(10-*''erg s-l) 


(lO^^ergs-i) 


(lO^^erg s-i) 


(keV) 


(keV) 




NSCS J121847+484410 


0.050 ± 0.006 


0.164 ± 0.018 


0.037 ± 0.020 




0.7 


SD 


NSCS J011502+002441 


0.404 ± 0.019 


0.508 ± 0.042 


0.229 ± 0.049 


r, c-o +0.09 

^•JJ -0.09 


2.1 


SD 


NSCS J100242+324218 


0.373 ± 0.023 


0.651 ± 0.040 


0.226 ± 0.032 


2.0 


SD 



Note. - A portion of this table is shown here for guidance regarding its form and content. A full version is available in the 
electronic edition of the MNRAS. 



X 

< 



(kcV) [ClliSl 

Figure 12. Ratio between the temperature measures available 
in CIRS and selected from the literature (BAX). The results are 
shown as a function of the CIRS values. The figure shows the 18 
common clusters with temperature available. 



from two sources only (IJones fc FormanI Il999l : iHorneij 
I2OOII ). while the latter comes from all the references listed 
above. As we see, the agreement for the 18 common clusters 
is very good, with almost all systems being concordant 
within 10%. 

For these cluster s, we employed the M200-TX relation 
given by equation 3 of lPopesso et al.l (|2005l ) to estimate mass 
from the temperature values. This relation was calibrated 
using a sample with more than one hundred clusters. Our 
goal is to compare the mass estimates obtained through dif- 
ferent, independent methods. In one case the mass is esti- 
mated from the virial analysis applied to the galaxy dis- 
tribution within clusters, while in the other case, mass is 
derived from the tight connection of the cluster mass to the 
temperature of the intra-cluster gas. Note that, from the 183 
low— 2; systems studied here, temperature is only available 
for generally massive clusters, with cr > 400 km/s or M200 > 
1O"M0. 



5 SUBSTRUCTURE ESTIMATES 

We have used photometric and spectroscopic data for galax- 
ies in the NoSOCS and CIRS clusters to obtain an esti- 
mate of the fraction of systems with evidence for substruc- 
tu re. The codes used fo r the substructure analysis are those 
of lPinknev et al.] l|l996l ). who evaluated the performance of 
thirty-one statistical tests. For the current work we consider 
the results of two of these tests, which are known to be the 
most sensitive. The first test considered is the DS or A test 
(jPressler &: ShectmanI [1988I ) . which is a three dimensional 
test. The algorithm computes the mean velocity and stan- 
dard deviation (cr) of each galaxy and its N„„ nearest neigh- 
bors, where N„„ = N^^^ and N is the number of galaxies in 
the cluster region. The mai n difference to the approach of 
iDressler fc ShectmanI (|l988l ) is that here we consider N„,i = 
N^/'^, instead of N„„ = 10. Then these local means and cr 
values are compared with the global mean and o (based on 
all galaxies). For every galaxy, a deviation from the global 
value is defined by the formula below. Substructure is esti- 
mated with the cumulative deviation A (defined by 5i; 
see below). For objects with no substructure we have A ~ 
N. 

^2 ^ ^ ^"^2"^ ^ ) [(""local - V)'^ + [criocal - Crf]- (8) 

The second test employed is two-dimensi onal, called 
symmetry or /3 test and it was introduced by IWest et al.l 
()1988. ). The test looks for significant deviations from mir- 
ror symmetry about the cluster center. It assumes that a 
substructure represents a local asymmetry superposed on 
an otherwise symmetric distribution. For every galaxy "i", 
a local density estimate di is obtained from th e mean dis- 
tance to the N^^'^ nearest neighbors (note that IWest et al.l 
Il988l used the five nearest neighbors). The local density, do, 
for a point "o" diametrical to a galaxy "i" is estimated in the 
same way. For a symmetric galaxy distribution the values of 
di and do should, on average, be approximately equal, but 
they will differ for clumpy distributions. The asymmetry for 
a given galaxy "i" is given by 



A^iog(| 



(9) 



The /3-statistic is then defined by the average value 
<Pi> over all galaxies. For a symmetric distribution < (5 >~ 
0, while values of </3> significantly greater than indicate 
asymmetries. 

It is well known that any substructure test statistic 
has little meaning if not properly normalized, which can 
be achieved by comparing the results for the input data to 
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those for substructure- free samples (the null hypothesis). For 
a two-dimensional test the null hypothesis is an azimuthally 
symetric smooth distribution, in which the surface density 
decreases with radius. For a three-dimensional test the null 
hypothesis is given by no correlation between position and 
velocity. In other words, velocity mean and dispersion should 
be the same locally and globally. For both tests, the signifi- 
cance level is determined through Monte Carlo simulations. 
For the /3 test the null hypothesis files are created through 
azimuthal randomizations of the positions. In that case, the 
distance of each galaxy to the cluster center is kept, while an 
azimuth is randomly assigned. For the A test the velocities 
are shufHed r andomly with respec t to the positions, which 
remain fixed (|Pinknev et al.lll99d ). 

For each input data set, we generate five-hundred sim- 
ulated realizations. We then calculate the number of Monte 
Carlo simulations which show more substructure than the 
real data. Finally, this number is divided by the number of 
Monte Carlo simulations. For the current work we set our 
significance threshold at 5%, meaning that only twenty-five 
simulated data sets can have substructure statistics higher 
than the observations to consider a substru cture estimate 
signifi cant. Further det a ils ca n be found in IPinknev et al.l 
(|l996t ) and iLopes et all (|2006l '). The values of the two sub- 
structure tests and their significance levels are listed in the 
last four columns of Tables f and 3 for the NoSOCS and 
CIRS samples, respectively. 

The substructure tests are only applied to clusters with 
at least five galaxies available within the aperture being 
used. Four different cases are considered. In the first we use 
the A test on the galaxies in the 170 clusters with at least 
five galaxies with available Zgpec within R2oo- We find that 
thirty-nine (~ 23%) clusters show signs of substructure. In 
the second case, we applied the /3 test to the same data set as 
above. This time, twenty-six (~ 15%) are found to have sub- 
structure. Third, we applied the /3 test to all galaxies within 
R200 in each cluster with no restriction for using only galax- 
ies that have Zspec available. Th e whole photometr ic data at 
m* — 1 ^ TTir ^m*-|-f is used (see lLopes et al.l2006l for details 
regarding this choice). For this case, we find that forty-three 
(~ 24%) clusters of all 179 show strong signs of substruc- 
ture. Finally, we applied the (3 test to all galaxies within 1.5 
h~^ Mpc of each cluster and found that sixty- four clusters 
(~ 35%), of all 183, show strong signs of substructur e. Thi s 
last result is in line with the findings of Lopes ct al. 
The rate of overlap between the four cases above is ~ 50%. 



6 COMPARISON OF OPTICAL AND X-RAY 
MASS ESTIMATES 

In this section, we compare mass estimates obtained from 
the dynamical analysis of the galaxy distribution with the 
M2oo-Tjf relation (Tx given in BAX). In Figures 13 and 14, 
we show the ratio between the optical and X-ray masses. 
First, we considere d the mass estimates o btained with the 
caustic technique l|Rines fc Diafericl |2006| ). That is shown 
in Figure 13 for the nineteen CIRS clusters with measured 
gas temperatures. Then we show in Figure 14 the results 
achieved when considering the mass estimates from this 
work. In both plots, we show clusters with substructure (es- 
timated with the A test) as open circles. The thin solid line 



on each figure shows the ratio equal to one, while the thick 
dashed line indicates the median value of the ratio for clus- 
ters with no substructure. 

The inspection of these figures shows more clusters in 
Figure 14 agreeing within 40% (the typical error in mass 
measurements) than in Figure 13. A few outliers displayed 
in Figure 14 (ratios lower than 0.5) are clusters affected by 
substructure. In this work, we find a good agreement with 
the X-ray expectations with a very weak trend in the sense 
of smaller X-ray to optical mass ratios for smaller masses. 
Moreover, we note that our result a grees to the one shown 
in Figure 1 of iPopesso et al.l (|2005h . where they also show 
a trend for optical estimates being slightly larger than the 
X-ray values. 

In Figure 13 the two clusters with mass ratios close to 
5 are Abell 1650 and MKW 08. We do not have a direct 
explanation for the lower value of M200 found by the caustic 
technique. After searching for other temperature measures 
in the literature, we found those to be consistent to the ones 
we obtained from BAX. We can also see in Figure 14 that 
the virial mass we determined is well matched to the X- 
ray based mass, for these two clusters. The virial and X-ray 
masses differ only by 10.4% for AbeU 1650 and 2.6% for 
MKW 08. The major outlier in Figure 14 is AbeU 2197. The 
inspection of its phase-space diagram (Figure 5) does not 
qualify it as a problematic cluster. We note that this is one 
of the clusters found to have substructure (by both tests, 
A and /3), which could explain the higher optical estimate. 
iRines fc Diaferid (I2OO6I ) comment that A2197 is composed of 
two X-ray groups (A2197W and A2197E). The X-ray center 
is considered the peak of A2197W and the optical centroid 
lie within the two X-ray groups. 

We also notice that in Figure 13 most of the clusters 
with substructure (4 out of 6) have mass differences within 
40%, while in Figure 14 they all have differences higher 
than 40%. From the comparison shown in Figure 14 we con- 
clude that clusters with substructure generally have larger 
mass differences. That conclusion is not supported by Fig- 
ure 13, indicating that the virial masses are more affected 
by projection or substructure, when compared to the caus- 
tic results. If we consider clusters with substructure to be 
not virialized the velocity dispersions generally overestimate 
the mass. The caustic mass is obtained directly from the 
mass profile, being independent of op. The absolute me- 
dian difference in Figure 13 is also higher than in Figure 14. 
However, the virial masses are generally higher than the X- 
ray masses. The opposite occurs when comparing X-ray and 
caustic masses. We also note there are more clusters closer 
to the ratio of one in Figure 14. 

Considering the results presented in Figures 13 and 14 
we carried out a mo re detailed comparison study with the 
results presented bv iRines fc Diaferid (|2006l ). As both esti- 
mates are obtained from the galaxy distribution, we would 
not expect large differences. However, note that the inter- 
loper removal procedure is very different. Here, we apply 
the shifting gapp e r tech nique bv lFadda et"al] (|l996l ). while 
iRines fc Diaferid l|2006l ) employ the caustic method. The 
maximum aperture defining the region where galaxie s are 
selected is also very distinct. iRines fc Diaferid (|2006l ) con- 
sider an aperture of 10 h~^ Mpc, four times the size con- 
sidered here. Finally, the most relevant difference in the 
two approaches is that our masses result from the applica- 
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Figure 13. The ratio between the mass obtained from the 
dynamical analysis o f the optical data (masses measured by 
iRines fc Diaferiol200d) and the mass estimated from Tx (selected 
from BAX). The figure shows the 19 common clusters with tem- 
perature available in BAX. Open circles represent clusters with 
substructure (estimated with the A test). The thin solid line rep- 
resents the ratio of one, while the thick dashed line shows the 
median value for clusters without substructure (filled circles). 



Figure 14. The ratio between the mass obtained from the dy- 
namical analysis of the optical data (masses from this work) and 
the mass estimated from Tx- The figure shows the 21 clusters 
with temperature available in BAX. Open circles represent clus- 
ters with substructure (estimated with the A test). The thin solid 
line represents the ratio of one, while the thick dashed line shows 
the median value for clusters without substructure (filled circles). 



tion of the virial theore m to the galaxy distribution, while 
iRines fc Diaferiol (|2006h estimate masses using the caustic 
technique (being R200 and M200 directly obtained from the 
caustic mass profile). 

In Figure 15, we show the comparison of the two optical 
mass estimates for the fifty-six clusters in common. Clearly 
many clusters disagree by more than 40%. Substructure can- 
not explain the higher masses obtained by the virial analysis 
as there are clusters exhibiting substructure and still agree 
well in mass. We also noticed that clusters with higher ratios 
(> 4) are generally the poorer ones and most of these have 
no indication of substructure. The higher masses obtained 
with the virial analysis are in agreement to what is seen in 
Figures 13 and 14, indicating a generally better agreement of 
our measurements to the X-ray estimates. The original CIRS 
values tend to be lower than the X-ray ones and thu s tend 
to be lower than our results. iRines fc Diaferiol l|2006l ) argue 
that they found a good correlation between M500 (caustic es- 
timate) and Tx , as well as between M200 (from the caustic) 
and Gp measured within R200 (note that their crp is inde- 
pendent of M200 and R200). They finally compare the caustic 
and virial masses at R200 finding an excellent agreement be- 
tween these measurements, thus concluding that the c aus- 
tic ma ss estimate is unbiased. Note that Rines fc Dia feriol 
(|2006l ) consider the R200 obtained by the caustic to deter- 
mine the virial mass. 

The disagreement we find with lRines fc Diaferiol (|2006l ) 
could result primarily from the interloper removal proce- 



dure and the velocity dispersion estimates. However, we see 
in Figure 16 that just a few systems have estimates of crp 
differing by more than 20%. Hence, we conclude that the in- 
terloper removal method does not affect many clusters. Note 
also that our ap values plotted are obtained only with clus- 
ter galaxies within R200. Those differ little from the ones ob- 
tained within the maximum aperture Ra (generally close to 
2.5 h~^ Mpc) and listed in Table 4. However, that does not 
significantly affect the estimates if Ra is sufficiently large (as 
seen in Figure 7). The good agreement exhibited in Figure 
16 leads us to conclude that the difference in mass estimates 
does not result from the interloper removal method. 

Hence, we compare in Figure 17 the values of R200. We 
see that our physical radius is generally different from the 
ones from the caustic. Namely, R200 determined from the 
caustic method is systematically smaller than our values. 
Thus, the assumptions we make in the virial analysis lead 
to intrinsically different radius and mass estimates than the 
ones obtained from the mass profile, derived from the caustic 
analysis. The main reasons for the discrepancies may actu- 
ally lie in the construction of the caustic mass profile, from 
which R200 and M200 are ulti mately derived. Details about 
the caustic profile are found in I Rines fc Diaferiol (|2006l ). One 
important conclusion these authors reach is that the caustic 
mass profile is better represented by an NFW profile for 50% 
of their sample. For 49% of the clusters a better fit is given 
by a Hernquist profile and for 1% by a singular isothermal 
sphere (SIS). In our approach we consider the NFW profile 
for representing the clusters. Depending of the concentration 
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Figure 15. Ratio between optical masses obtained by 
iRines fc Diaferiol ||2006| ') and in this work. The figure shows the 56 
CIRS clusters used here. The thin solid line represents the ratio 
of one, while the thick dashed line shows the median value for 
clusters without substructure (filled circles). 



parameter the NFW profile may be lower than the caustic 
mass profile at low radius, resulting in a smaller R200 value 
estimated by the c austic (see their F ig ure 12 ). 

In summary, iRines fc Diaferiol l|2006l ) derive mass 
from the caustic mass profiles and they claim concor- 
dance between the caustic results and their virial masses. 
However, comparing o ur masses and the caustic values 
mes fc Diaferiol I2OO6D to the X -ray estimates, we find a 
better correlation with our virial masses. Therefore, we con- 
clude that our estimates better represent the cluster poten- 
tials. That conclusion is also corroborated by the scaling 
relations we show in paper II (|Lopes et al.ll200a ). where we 
find more accurate results when usi ng our mass estimates 
than the ones fi-omlRines fc Diaferiol (2006). It is important 
to stress that iRines fc Diaferiol (|2006i ) use R200 determined 
by the caustic to compute the virial masses. That choice 
contributes for the good agreement they find between the 
caustic and virial masses. 



Figure 16. Same as previous figure, but showing the results re- 
garding the velocity dispersions. 



7 CONCLUSIONS 



We have used SDSS data to study the galaxy properties 
of clusters first selected in the optical using DPOSS. We 
extr acted SDSS data for every NoS OCS supplemental clus- 
ter (|Lopedl200j : iLopes et al]|2004 ) sampled in DR5. Using 
the high quality photometric data o f SDSS, we estimated 



new photometric redshi fts (following LopesllioOTi ). richness 
and optical luminosity ijLopes et al.ll2006l ). After removing 
potential double clusters, we compiled a list with 7,414 sys- 
tems (~ 75% of the original NoSOCS supplemental sample) 
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Figure 17. Same as previous figure, but showing the results re- 
garding the physical radius (R200). 



with properties derived from SDSS data. This cluster list 
has systems to 2 ~ 0.5. When estimating richness, we also 
assess the effect of — correction and evolution in M* to 
the richness estimates. 

For a subset of systems a.t z ^ 0.10 we made use of the 
SDSS spectroscopic data to identify groups in redshift space. 
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confirming tliese cluster candidates as true 3D overdensi- 
ties. Fo r every confirmed s ystem, we applied the "shifting 
gapper iFadda et al ] (ll99d) technique to exclude galaxy in- 
terlopers in the phase-space distributions and estimate the 
velocity dispersions and masses (M500 and M2oo)- The fi- 
nal sample with these properties derived at low-z comprises 
127 clusters. As the NoSOCS supplement does not sam- 
ple very massive clusters at low— 2, we extended the anal- 
ysis for richer systems including those clusters from CIRS 
iRines fc Diaferioll2006l ). Optical and dynamical properties 
are derived for fifty-six CIRS systems. Hence, the final list 
of clusters studied at low redshift (2 ^ 0.10) comprises 183 
objects (127 from NoSOCS and 56 from CIRS). The main 
results of this work are: 

(i) We have checked how velocity dispersion and mass are 
affected by a less rigorous method for interloper removal, 
as well as the radial cutoff (maximum aperture sampled) 
and the luminosity limit of the spectroscopic survey. We 
found that looser criteria for removing interlopers can lead 
to higher mass estimates, which could bias the mass cali- 
bration (see discussion in paper II). With regards the max- 
imum radius, we found that an aperture of 2.5 h~^ Mpc is 
large enough to give unbiased velocity dispersion estimates. 
However, we notice that op and mass estimates are severely 
biased if one considers a spectroscopic survey which is not 
complete to at least M* -I- 1. That bias affects the construc- 
tion of cluster s caling relations. 

(u) Following ^al et al.l l|2008l ) we used RASS data to es- 
timate X-ray luminosities for all 183 clusters. We obtained 
reliable X-ray detections for 85% of the clusters, while for 
the remaining only upper limits were measured. 

(iii) For a subset of these clusters (21 objects), we se- 
lected temperature measurements from the literature and 
estimated cluster mass assuming an M-Tx relation. The 
comparison of optical and X-ray mass estimates for the sub- 
set of twenty-one clusters with measured Tx reveals good 
concordance between masses from different wavelengths. We 
note a small trend for the optical masses estimated in this 
work to be larger than the X-ray values. However, we found 
the cluster masses dete rmined from the caustic technique 
IIRines k, Diaferid[200^ ) show a more pronounced trend for 
lower masses compared to masses determined from the gas 
temperatures. Thus, the optical masses from the caustic 
technique are generally smaller than the masses from the 
virial analysis obtained in the current work (which we con- 
clude better represents the cluster potential). 

(iv) When comparin g the velocity di s persio n estimates 
of the present work to iRines fc Diaferiol (|2006l ) we do not 
see such a difference, which indicates that the interloper re- 
moval procedures lead to similar lists of cluster members. 
The comparison of R200 values determined in this work and 
by iRines & Diafcrio (2006) reveals a systematic difference, 
which may explain the difference in mass. The mass discrep- 
ancies may be due to the shape of the caustic mass profile, 
which may result i n biased values of R200 and thus M200 . In 
iLopes et al.l (|2008l ) we have also shown that the virial masses 
we obtained lead to more accurate scaling relations than the 
ones based on the caustic approach. More important, the 
scaling relations found with X-ray masses fully agree with 
the virial masses based relations, but show distinct results 
to the relations derived from the caustic masses. The scaling 



relations based in the caustic masses also show a much larger 
scatter than thos e based in the virial masses (see paper II; 
iLopes et al]|2008l ). 

(v) For all 183 systems we have also estimated substruc- 
ture from the galaxy distribution in two, and three dimen- 
sions, using the /3 and A tests, respectively. We investigated 
the effect of substructure on the mass estimates. When con- 
sidering only the spectroscopically confirmed cluster mem- 
bers, we find that 23% of clusters have substructure esti- 
mated from the 3D A test, while 15% are found to have 
substructure with the 2D j3 test. When considering the pho- 
tometric data set we find that 24% of clusters have sub- 
structure according to the /3 test within R2oo- The same 
test applied to the photometric data within 1.5 h~^ Mpc 
estimates in 35% the ratio of systems with substructure. We 
also find that most of the clusters with substructure show 
large differences in mass estimates, when comparing X-ray 
and virial (this work) measures. The same is not true if the 
X-ray masses are compared to the caustic values. We con- 
clude that substructure can partly explain the differences in 
the mass values from optical and X-ray data. 
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